The role of inter-well tunneling strength on coherence dynamics of two-species 

Bose-Einstein condensates 
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Coherence dynamics of two-species Bose-Einstein condensates in double wells is investigated in 
mean field approximation. We show that the system can exhibit decoherence phenomena even 
without the condensate-environment coupling and the variation tendency of the degree of coherence 
depends on not only the parameters of the system but also the initial states. We also investigate the 
time evolution of the degree of coherence for a Rosen-Zener form of tunneling strength, and propose 
a method to get a condensate system with certain degree of coherence through a time-dependent 
tunneling strength. 
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I. INTRODUCTION 

The manipulating Bose-Einstein condensates (BECs) 
in double wells provides us a versatile tool to explore the 
underlying physics in various nonlinear phenomena since 
almost each parameter can be tuned experimentally 
There have been many studies on the fascinatingfea- 
tures of nonlinear effect, such as Rabi oscillation 043; 
Josephson oscillation (^-{j, self trapping d, [7l-[l0|. and 
measure synchronization 11], in terms of BECs in dou- 
ble wells. Because the BECs in double wells can be re- 
garded as a two-level system, it is also expected either 
to be employed as a possible qubit or to simulate certain 
issues [12l] in quantum computation and information. Re- 
cently, the effect of decoherence of BECs in double wells 
was investigated experimentally by means of interference 
between BECs [l^, |l3| and studied theoretically in terms 
of single-particle density matrix [l5| . In order to ex- 
hibit the phenomena of decoherence, they [T3l - [l5j need to 
introduce the condensate-environment coupling because 
one-species BEC in double wells were merely considered 
there. In comparison to one-species BEC system, the 
two-species system can exhibit distinct decoherence phe- 
nomena due to the existence of the interspecies inter- 
action. For example, the degree of coherence in a two- 
species system can evolve with time without the applica- 
tion of condensate-environment coupling, which is useful 
for one to get a system with desired degree of coherence. 
It is therefore worthwhile to study the coherence dynam- 
ics of two-species BECs in double wells. 

In this paper, we study the effects of the inter-well 
tunneling strength on the coherence dynamics for a sys- 
tem of two-species BECs in double wells. With the help 
of the reduced single-particle density matrix, we show 
that such a system can exhibit decoherence phenomena 
without condensate-environment coupling. We also pro- 
pose an experimental strategy to prepare a BEC system 
with any desired degree of coherence through a time- 
dependent tunneling strength. In the next section, we 
model the two-species BEC system and introduce the re- 
duced single-particle density matrix. In Sec. lIIIl we study 



the time evolution of the degree of coherence for a time- 
independent inter- well tunneling strength. In Sec. lIVl we 
investigate the time evolution of the degree of coherence 
for a Rosen-Zener form of tunneling strength and discuss 
the possibility of preparing a BEC system with any de- 
gree of coherence. Then we briefly give our conclusion in 
Sec. El 



II. MODEL AND REDUCED DENSITY 
MATRIX 

We consider a two-species Bose-Einstein condensate 
system confined in double wells. The Hamiltonian is 
given by 



H = ~Ja{a\a2+ll.c.)~ Jbib{b2+R.C.) + J2Ua 



b^ai^bi 



-zUaaflaiflai + ^ -zUbbfl-bi 



n-bi 



(1) 



where a| (6|) and (hi) creates and annihilates a bosonic 
atom of species a (6) in the ith well, respectively; hai — 
alcii {fibi = b\bi) denotes the particle number operator 
of species a (b). Here the parameters Ja and Jb de- 
note the tunneling strengths of species a and b between 
the two wells, Uaa and Ubb are the intraspecies inter- 
action strengths, and Uab is the interspecies interaction 
strength. The Hamiltonian ([IJ can describe a BEC mix- 
ture confined in a double well potential consisting of dif- 
ferent atoms, or different isotopes, or different hypcrfine 
states of the same kind of atom. The coherence dynamics 
of the above model has not been investigated although its 
dynamical properties, like Josephson oscillation, stability 
and measure synchronization etc., have been studied in 
earlier works [1, [ll|, [l^. Whereas, we know that the 
role of coherence of the system is very important since 
the first obstacle attempted to be avoided is the decoher- 
ence when a condensate in double wells is expected to be 
employed as a qubit. So the coherence dynamics of the 
model (HI) is worthy of study. 
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As we know, under the semiclassical limit 0, [13] , 
the dynamics of this system is conventionaUy studied in 
mean-field approach by replacing the expectation values 
of annihilators with complex numbers, i.e., { ai ) — hi 
and { bi) = bi. With the help of Heisenberg equation of 
motion for operators, one can easily get the dynamical 
equations for a,; and bi. These equations guarantee the 
conservation law |aip + |a2p — N/2 and |6ip + |62p = 
N/2 with N being the total particle number of species a 



and b. To simplify the calculation, one usually assumes 
Ui — (ii/y/N and bi = bi/^/N. Then the aforementioned 
dynamical equations for di and bi can be rewritten as, 

^J^m^H,ff\^P) (2) 

with 



Hcff 



Ua\ai\^ + Uab\bl\' 
-Ja 






-Ja 

Ua\a2\^+Uab\b2\^ 






Ub\bi\' 






+ Uab\ai\ 



Ub\b2 







-Jb 

+ Uab\a2\^ J 



(3) 



where Ua = UaaN, Ub = UbbN , and Uab = UabN. Here 
the wave function | ^) refers to 

li/') = ( ai, 02, fei, h2)^ = Y^ \(7)® (4) 

(T.i 

where the state | cr) =| a) or | &) specifies the two different 
species while | i) =| 1) or | 2) specifies the two wells. 
Note that the dynamical properties of the system can 
be determined by Eq. ^ if the system is in a completely 
coherent state (i.e., a pure state). Whereas, if the system 
is in a mixed state, the equation ([2]) becomes insufficient. 

In order to study the coherence dynamics of the sys- 
tem, we introduce the single-particle density matrix p = | 
'>p){4' I whose elements are p^^(^, = 1,2,3,4). From 
this definition, we can see that, as a 4 x 4 matrix, p de- 
scribes a pure state. Their diagonal elements pn (P22) 
and (P44) represent the population of species a (b) 
in the first and second well, respectively. In this paper, 
we only focus on the distribution of the total particle 
numbers in the two wells but not distinguish the particle 
species, so the system can be described by the reduced 
density matrix 

p^ = J2{a\p\a), (5) 

cr 

whose elements are p^^ = pu + P33, p\2 = P12 + P34, 
P21 = P21 +P43, and P22 — P22 + P44- Clearly, the matrix 
fF can describe a mixed state. Its diagonal elements pii 
and P22 represent the total population probability in the 
first and second well, respectively. To investigate the 
coherence dynamics of the system, we can introduce the 
definition of degree of coherence according to Ref . d, |T3] , 

^ = Tripn' - 1. (6) 

From Eq. ([SJ and Eq. ©, we can see that the time evo- 
lution of the degree of coherence depends on that of the 



elements p^i/ of the single-particle density matrix p. Since 
p describes a pure state, the time evolution of can 
be determined by Eq. (0). Thus one can solve Eq. ^ 
to get the evolution of Ui and bi firstly, and then gives 
the time evolution of p^j/ according to the definition of 
p — \^){'ip\. Then in the following sections, we will study 
the coherence dynamics of the system with the help of 
Eq. (2)). Note that in the following calculations, we take 
Ja = Jb = J for simplicity without losing the generality. 



III. EVOLUTION OF DEGREE OF 
COHERENCE 

We know that the degree of coherence does not change 
for an isolated BEC system Q. To study the effect of 
decoherence of BEC system, several authors considered 
the condensate-environment coupling [l^ll^. Whereas, 
we will show that the degree of coherence of the two- 
species BEC system in double wells can still change with 
time even without the condensate-environment coupling. 
Here we consider the case of time-independent inter- well 
tunneling strength, i.e., J is a constant in the calcula- 
tion. Due to the fact that Eq. ([2]) can not be analytically 
solved, we solve Eq. ^ numerically to get the time evo- 
lution of Qi and bi. Then we give the time evolution of 
the degree of coherence with the help of the definition 
of p and pr and Eq. The corresponding results are 
summarized in Fig. [TJ 

In Fig. [TJ we plot the time evolution of the degree 
of coherence for different initial states and parameters. 
The initial states are ai — \J\/2, 02 = 0, &i = ■\/l/2 
and 62 = for the blue (top) line, a\ = ■\/2/5, 02 = 
yr/lO, bx = ^/Tj2 and 62 = for the red (middle) line, 
and ai — \J\/2, 02 = 0, 61 = and &2 = \/l/2 for the 
black (bottom) line. From the Fig.[T](a), we can find that 
the degree of coherence of the two-species BEC system in 
double wells does not change for any initial states when 
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FIG. 1: (Color online) Time evolution of the degree of co- 
herence for different initial states. The parameters are Ua = 
Ub = Uab = 3 J (a), Ua = Ub = 3 J and Uat = 6 J (b), 
Ua = 3 J, Ub = 2 J and (/at = J (c), and Ua = 3 J, Ub = 2 J 
and [/„(, = 6 J (d). 



and the parameters of the system, so one can change it 
by varying either the initial states or the parameters of 
the system. In the following, we will show how to control 
the degree of coherence through a Rosen-Zener form of 
J, 



J^jQs'm^ujt (0 < i < tt/w), 



(7) 



i.e., J increases from zero to its maximum value Jq and 
then decreases to zero again in the end of the calculation. 
In the numerical calculation in this section, we take the 
initial state ai — VOA, 02 = VO.l, bi — \/0.5, and 62 = 0. 




Ua — Ub — Uab- That can be easily understood because 
this two-species system is equivalent to the one-species 
system once Ua = Ut = Uab- And according to Ref. ([l]), 
the degree of coherence of the one-species BEG system 
in double wells does not change without the condensate- 
environment coupling. So the result shown in Fig. [T] (a) 
is reasonable. The figure [T] (b) shows that when Ua = 
Ub ^ Uab, the degree of coherence does not change for 
the case of a; — hi at the initial time, i.e., the initial 
population distributions of species a and b are the same, 
but that changes for the other initial states. Note that 
for the case oi Ua = Ub 7^ Uab, the time evolutions of 
species a and b are the same if = bi at the initial 
time. So the degree of coherence does not change due to 
the two species having the same symmetry. Additionally, 
from the Fig. [1] (c) and (d), we can find that the degree 
of coherence changes for any initial states for the case of 
Ua ^ Ub Uab- Comparing Fig.[T](c) and (d), we can see 
that for the same initial state, the variation tendency of 
the degree of coherence depends on the parameters of the 
system. In summary, whether the degree of coherence 
can change with time depends on both the parameters 
and the initial states. So does the variation tendency of 
the degree of coherence. This fact indicates that one can 
control the degree of coherence without introducing the 
effect of environment, which will be discussed in the next 
section. 



IV. CONTROLLING THE DEGREE OF 
COHERENCE 

In previous section, we show that the time evolution 
of the degree of coherence depends on the initial states 



FIG. 2: (Color online) Time evolution of the degree of coher- 
ence for the Rosen-Zener form of inter- well tunneling strength. 

2 11 
The parameters are Ub = -Ua,Uab = -^Ua,Jo = "^UayUj = 

O.OlUa. 
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FIG. 3: The dependence of the maximum value of the degree 

of coherence (left panel) and the corresponding time (right 

panel) on the period of the inter- well tunneling strength. The 

2 11 
parameters are Ub = —Ua, Uab = —Ua, Jo = -Ua. 



The time evolution of the degree of coherence for the 
Rosen-Zener form of inter- well tunneling strength is plot- 
ted in Fig. O We can see that the degree of coherence 
changes with time and can reach the maximum value 
??max at time t = tmax- Notc that the values of yymax and 
t = imax depend on both the period w and the maximum 
value Jo of the inter-well tunneling strength. In Fig. |31 



we plot the dependence of ?]„ 
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FIG. 4: The dependence of the maximum value of the de- 
gree of coherence (left panel) and the corresponding time 
(right panel) on the maximum value of the inter-well tunnel- 

2 1 

ing strength. The parameters are Ub = g'^"' ~ g'^i, Jo = 
it/.. 



Jo in the left and right panel, respectively. Meanwhile, 
we plot the dependence of r^max and t,„ax on Jq for a fixed 
u! in Fig. 01 

From Fig. [3] and Fig. |4l we can see that the values 
of Lo and Jq affect the maximum value of the degree of 
coherence sensitively. In order to get a system with large 
value of degree of coherence, one must choose the form of 
J with suitable period uj and maximum value Jq. Note 
that although one can get a system with large value of 
degree of coherence through tuning the value of w and Jq , 
it is difficult to control the time evolution of the degree 
of coherence due to the fact that the degree of coherence 
does not evolve periodically, which can be confirmed by 
Fig. [2] In order to overcome the aforementioned problem, 
we consider the following form of J 



J 




(0 < t < imax) 



(8) 



where i,nax can be obtained from Fig. [3] and Fig. |31 

For the form of J given in Eq. ([S]), the time evo- 
lution of the degree of coherence is plotted in Fig. [S] 
From this figure, we can find that the degree of coher- 
ence reaches its maximum value which is determined 
by the values of uj and Jq, and then oscillates period- 
ically at the following time. Since Eq. ^ can be an- 
alytically solved for the case of J = 0, we can give 
the oscillation period in analytical method. Solving 
Eq. ([2]), we obtain ai{t) — ai(tniax) exp (— iA^i), and 
bi{t) = 6i(tmax)exp(-jBit) where Ai = ?7o|ai(imax)P + 
Uab\h{t„,a^)\^ , and B, = Utlhit^^^)]"^ + C/af,|aj(imax)P- 
Substituting the expressions of Ui and bi into the Eq. ([SJ , 
we can get the oscillation period of the degree of coher- 
ence T = 2tt/\A2 — Ai + Bi — B2\. For the parameters 
taken in Fig. [SJ we find TUa — 56.5 which is consistent 
with the numerical result. From the above discussion. 




FIG. 5: (Color online) Time evolution of the degree of co- 
herence for the inter- well tunneling strength given in Eq. 

2 11 
The parameters are Ut = -Ua,Uab = -^Ua,Jo = -^Ua,ijJ = 

0.005t/a. 



we know that the degree of coherence oscillates period- 
ically after time tmax, and both the maximum value of 
the degree of coherence and the oscillation period can 
be changed by tuning the values of oj and Jq. Once the 
values of lo and Jq are fixed, the time evolution of the 
degree of coherence is well defined, so that we can know 
the degree of coherence of the system at any time. Then 
one can easily get a system with any desired degree of 
coherence through controlling the evolution time. 



V. SUMMARY 

In the above, we investigated the coherence dynam- 
ics for a system of two-species Bosc-Einstein conden- 
sate in double wells. In mean field approximation, we 
studied the infiuence of the inter- well tunneling strength 
on the coherence features with the help of the reduced 
single-particle density matrix. Since we need not dis- 
tinguish particle species in the system, we only focused 
on the distribution of the total particle numbers in the 
two wells, which can be described by a 2x2 reduced den- 
sity matrix. After studying the time evolution of the de- 
gree of coherence for a time-independent inter-well tun- 
neling strength, we found that the degree of coherence 
of the two-species BEG system changes for some pa- 
rameters and initial states even without the condensate- 
environment coupling, which differs from the case of 
Refs. [H, Motivated by the fact that the variation 
tendency of the degree of coherence depends on both the 
parameters and the initial states of the system, we con- 
sidered a system with a Rosen-Zener form of inter-well 
tunneling strength to control the degree of coherence. Al- 
though its tendency is not periodical, the degree of coher- 
ence can reach a maximum value 77max at time <max which 
are dependent of the period uj and the maximum value 
Jo of the inter-well tunneling strength. The dependence 
of ?7max and imax On UJ and Jo we obtained is helpful for 
one to get a system with large value of degree of coher- 
ence utilizing a Rosen-Zener form of inter-well tunneling 
strength. We also gave a useful form of the inter- well tun- 
neling strength for one to easily get a system with any 
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degree of coherence by controlling the evolution time. 10874149 and No. 11074216. 
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